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Abstract 

The metal-insulator transition observed in the In/Si(lll)-4x 1 reconstruction is studied by means of ab initio 
calculations of a simplified model of the surface. Different surface bands are identified and classified according to 
their origin and their response to several structural distortions. We support the, recently proposed [New J. of Phys. 
7 (2005) 100], combination of a shear and a Peierls distortions as the origin of the metal-insulator transition. Our 
results also seem to favor an electronic driving force for the transition. 



Quasi one-dimensional reconstructions formed 
by metal deposition on Si(lll) have been in- 
tensively studied in recent years. The electronic 
correlations and the coupling between electronic 
and structural degrees of freedom are enhanced 
in one-dimension and, as a consequence, several 
electronic and structural phase transitons are 
observed in these systems as the temperature 
is decreased. A nice example of this behavior is 
found in the In/Si(lll) system, which exhibits a 
4x1 — » 4x2 — » 8x2 structural transition accom- 
panied by a metal-insulator electronic transition. 

The room-temperature (RT) 4x1 structure of 
the In/Si(lll) surface is well established [1,2]. It 
consists of two neighboring zigzag In wires along 
the [112] direction. Each wire contains two In 
atoms per 4x1 cell and each In atom is bonded to 
one Si atom of the substrate. This model has been 
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confirmed by ab-initio calculations [3,4,5] which 
reproduce the scanning tunneling microscopy 
(STM) images [6,7], and the main features of the 
band structure. 

At RT the system presents three metallic sur- 
face bands with similar dispersion [8]. However, 
when the temperature is lowered below ~ 130 K [9] 
photoemission shows the formation of a band gap. 
This transition is accompanied with a doubling 
of the unit cell in the STM images [9]. The low 
temperature (LT) phase has been widely stud- 
ied experimentally [10,11,12,13]. Although several 
results seem to favor a model where the indium 
wires suffer a strong dimerization and break-up 
into trimers, the question is still far from settled. 
In fact, most ab-initio calculations find metallic 
structures which are only slightly distorted with 
respect to the RT phase [5,14]. Thus, the driving 
force behind the combined metal-insulator and 
structural transition remains unclear. 

Recently, an interesting mechanism for the 
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gap opening has been proposed by Ahn et al. in 
Ref. [15]. The occupation of the surface bands in 
the 4x1 structure is quite close to two electrons. If 
one of these bands is depopulated (the upper one), 
the other two become very close to half-filled and 
thus are suitable to suffer a Peierls transition due 
to a periodicity doubling. If this is true, it seems 
to indicate that: (i) there are at least two types 
of surface bands that originate in different regions 
of the substrate or have different symmetries, and 
(it) the metal-insulator transition is the result of a 
combination of two distinct structural distortions 
that couple with different bands. This last point 
is consistent with the recent first-principles cal- 
culations by Gonzalez et. al. [16]. These authors 
find an insulating 4x2 structure, reminiscent of 
that proposed by Kumpf et al. [1] , as a result of a 
combined shear and Peierls distortion of the 4x1 
RT phase. 

The objective of the present work is to under- 
stand in detail the origin and characteristics of the 
different electronic states involved in the metal- 
insulator transition, and how they couple to dif- 
ferent structural distortions. The emphasis is on 
the electronic bands associated with the indium 
atoms in the substrate. We use a simplified model 
(shown in Fig. 1) that captures the essence of the 
system. Our results support the main conclusions 
of Ref. [16], and point to a primary electronic ori- 
gin of the structural transition. Of course, in the 
real surface we can expect a delicate competition 
between the gain of electronic energy and the elas- 
tic energy associated with the different distortions. 
Work along this line is in progress and will be pub- 
lished elsewhere [17]. 

Our calculations were performed with the 
SIESTA method [18]. We used the generalized gra- 
dient approximation [19] to the density functional 
theory, norm-conserving pscudopotentials [20], 
and a double-^ polarized basis set of numerical 
atomic orbitals with cutoff radii corresponding to 
an energy shift of 200 mcV [18]. 1 The density of 
the auxiliary real-space grid [18] was equivalent to 
a plane- wave cutoff of 100 Ry. We typically used 12 
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inequivalent k-points along the wire direction (up 
to 250 in convergence tests, and 100 in the calcula- 
tions shown in Fig. 4). We used periodic boundary 
conditions with a vacuum distance of 15 A. 

Figure 1 presents the simple model used here 
to study the electronic properties of indium wires 
in the In/Si(lll) surface. We keep the two zigzag 
In wires present in the 4x1 unit cell, and substi- 
tute the neighboring silicon atoms with hydrogen. 
The In-H distances were optimized (1.86 A), and 
the relative positions of the In-H pairs were kept 
fixed for the rest of the study. This system pro- 
vides a qualitative model that, however, retains the 
main features of the electronic band structure of 
the In/Si(lll) surface (see below). The In-H bond 
is slightly more ionic than the In-Si bond. How- 
ever, we have checked that saturating with SiH3 
groups instead of H atoms leads to a very similar 
band structure. 

Figure 2 shows the evolution of the band struc- 
ture of our model as a function of the wire-wire 
distance d. There are two types of indium atoms 
in each wire, In(l) and In(2). The coordination of 
the In(2) atoms changes with d. Panel (a) corre- 
sponds to non-interacting (large d) zigzag wires, 
while panel (c) (rf=2.15 A) corresponds to a config- 
uration similar to that found in the In/Si(lll)-4x 1 
reconstruction. Although not completely evident 
due to appearance of interaction gaps, the band 
structure in Fig. 2 (a) can be rationalized in terms 
of three bands: (I) a strongly dispersive band as- 
sociated with the In(l) atoms and the In(2)-p y or- 
bitals (circles) , a flat band (II) with a clear contri- 
bution from In(l) and the In(2)-p 2 orbitals (trian- 
gles), and another flat band (III) with large In(2)- 
p x character (squares) . Taken into account the hy- 
drogen saturation, each In atom contributes with 
two valence electrons. Thus we have four electrons 
two distribute in these bands. Band (II) is doubly 
occupied and does not play any role in the argu- 
mentation below. Bands (I) and (III), however, are 
half-filled. Band (III) can be associated with the 
"dangling-bonds" in the In(2) atoms that project 
approximately into the x-direction and the vac- 
uum. 

As d is reduced, the interaction between the 
wires modifies the band structure. This can be seen 
in Fig. 2(b) and (c). Particularly, the dispersion of 
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bands derived from band (III) largely increase as 
a result of the overlap of the dangling-bonds in the 
neighboring wires. Finally, the electronic states as- 
sociated with this band become highly delocalized 
in the region between the two zigzag wires, and the 
band exhibits an almost free-electron dispersion. 
In the following we call this band the "interaction 
band" . 

Figure 2 (d) presents the band structure of the 
In/Si(lll)-4x 1 reconstruction calculated with 
a slab containing four silicon bilayers [17]. Four 
surface bands can be located in the gap of the 
silicon substrate. Three of them cross the Fermi 
level with similar dispersions. This is qualitatively 
reproduced by our model. The agreement is im- 
proved if a small shift of the Fermi level is allowed 
(see Fig. 2 (c)). The information from our model 
allows to catalogue the surface bands in two types: 
(i) two bands that are derived from the band (I) 
of the right and left wires, and (ii) the interaction 
band. The first two bands have a larger weight 
inside each of the wires and are quite sensitive to 
the structure of the zigzag chains. In contrast, the 
interaction band is localized in the region between 
the two zigzag wires and thus is more influenced 
by the relative positions of the wires. 

This division allows to envision a two-step 
route for the observed metal-insulator transition. 
Step one: the dimerization of the dangling-bonds 
from neighboring zigzag chains opens a gap in 
the interaction-band. This effect can be obtained 
without doubling the periodicity of the system, 
the relative displacement of the wires along their 
axes (see Fig. 3) suffices. Notice that this corre- 
sponds to the shear distortion in Ref. [16]. Step 
two: the remaining metallic bands need to ac- 
commodate two valence electrons. Since these two 
bands have very similar dispersions, both become 
approximately half-filled. As a consequence, the 
system is now suitable to suffer a Peierls transi- 
tion. This two-step mechanism is consistent with 
the calculations of Ref. [16] and the experimental 
evidence in Ref. [15]. 

Figure 3 (a) shows the evolution of the energy 
and the band structure of our system as function 
of the shear distortion. One of the zigzag wires 
was displaced along the y-direction by a magni- 
tude Ay with respect to the other. For each dis- 



placement the distance between the wires d was 
optimized. A gap is opened in the interaction band 
which widens with increasing Ay. For distortions 
larger than Ay ~0.5 A the Fermi level enters in this 
gap. This is reflected in the behavior of the energy 
that starts to decrease at this point. This behav- 
ior translates in an energy barrier of ~5 meV per 
In atom. The system is then left with two metallic 
bands that cross the Fermi level at nearby points 
in reciprocal space. This is also the case for the real 
In/Si(lll)-4x 1 surface as can be seen in Fig. 3 (b). 
However, in this case the behavior of the total en- 
ergy is different. Although still in the range of a 
few meV per In atom for moderate distortions, the 
shear deformation always increases the energy of 
the system [17]. 

We now study the effect of doubling the period- 
icity along the wires. We consider a quite simple 
Peierls-like distortion: the length of one every four 
bonds is shorten (the undistorted bond length is 
3.045 A). The distortion is identical for both wires. 
Still we have four different possibilities according 
to the different relative locations of the distorted 
bonds in the neighboring wires. This is illustrated 
in Fig. 4 (a). While doubling the unit cell, the four 
distortions break the symmetry of the system in 
different ways. This detail is quite important. The 
distortions open a gap at ^ due to the periodicity 
doubling. However, since the Fermi points do not 
exactly lie in that position, this does not guaran- 
tee that the system will become semiconducting. 
This is more clear for the extreme shear distortion 
(Ay=2.15 A, a=90°). In this case there is a mir- 
ror plane parallel to the axis of the wires which 
is only preserved by "distortion 1". As a conse- 
quence of this symmetry, the band structure of the 
system submitted to the "distortion 1" presents 
a band crossing and the system remains metallic. 
For the other three distortions the band crossing 
is avoided due to the break of symmetry and a 
gap opens at the Fermi level. This can be appreci- 
ated in Fig. 4 (b). For a ^90° the symmetry gap 
is always opened, although its magnitude depends 
again on each particular structural distortion. Fig- 
ures 4 (b) and (c) present the total energy for the 
different Peierls-like distortions (a shear deforma- 
tion has been previously applied to the system). 
"Distortion 2" and "distortion 4" are always the 



3 



most favorables. It is worth noting that the struc- 
ture proposed in Ref. [16] can be understood as the 
result of applying a combination of a shear distor- 
tion and the Peierls-like "distortion 2" presented 
above. 

In summary, we have studied the electronic 
structure of the indium zigzag wires seen on the 
In/Si(lll)-4xl reconstruction. The different sur- 
face bands are identified and classified according 
to their origin and their response to different struc- 
tural distortions. We confirm that the combination 
of a shear and a Peierls distortion, proposed in 
Ref. [16], provides a reasonable and robust route 
for the observed metal-insulator transition in this 
system. Our results also point to an electronic 
driving force of this transition. 
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Fig. 1. Our model of the In wires in the In/Si(lll)-4xl 
surface: it contains two zigzag indium wires saturated with 
hydrogen. 



Fig. 2. (a)-(c) Band structure (along the wires axes) of 
the system shown in Fig. 1 as a function of the wire-wire 
distance d: (a) isolated zigzag In wire (d=10.8 A), (b) 
d=4.31 A, and (c) d=2.15 A. Different symbols indicate the 
distinct character of the bands as determined from a Mul- 
liken population analysis [21]: In(l) and In(2)-p H (circles), 
In(2)-pa; (squares), and In(l) and In(2)-p z (triangles), (d) 
Band structure of the In/Si(lll)-4x 1 reconstruction along 
the [112] direction calculated with a slab containing four 
silicon bilayers [17], the bands with strong indium char- 
acter arc highlighted with circles. Energies are always re- 
ferred to the Fermi level. 



Fig. 3. Shear distortion, (a) Energy per indium atom 
as a function of the relative displacement Ay of the in- 
dium wires along their axes (see the scheme in the upper 
part of panel (a)). The band structures for three different 
values of Ay are also shown, (b) Band structure of the 
In/Si(lll)-4xl reconstruction with Ay=1.65 A. The inset 
shows the Brillouin zone, TX and YM run along the In 
wires. Energies in the band structures are referred to the 
Fermi level. 



Fig. 4. Peierls distortion, (a) Scheme showing four inequiv- 
alent Peierls-like distortions: the lenght of one of the bonds 
(indicated by an arrow) is modified by the same amount 
in both wires; different distortions correspond to different 
relative positions of the distorted bonds and are numbered 
according to the labels of the different bonds in the right 
wire. Panels (b) and (c) show the total energy per indium 
atom as a function of the modified bond length for dis- 
tortions 1 (solid), 2 (dashed), 3 (dotted), 4 (dash-dotted). 
Panel (b) corresponds to a=90° (Ay=2.15 A), while in 
panel (c) a=154° (Aj/=1.05 A). The insets show the band 
structures close to the Brillouin-zone boundary for distor- 
tions 1 and 2 in panel (b), and 2 in panel (c). Energies in 
the band structures are referred to the Fermi level. 
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